Strong Semiclassical Approximation of Wigner Functions for the 

Hartree Dynamics 

A. Athanassoulifl T. Paui, F. Pezzotti, M. Pulvirentfl 

Abstract 

We consider the Wigner equation eorresponding to a nonlinear Schrodinger evolution of the Hartree type 
in the semiclassical limit ft — ^ 0. 

Under appropriate assumptions on the initial data and the interaction potential, we show that the Wigner 
function is close in to its weak limit, the solution of the corresponding Vlasov equation. The strong ap- 
proximation allows the construction of semiclassical operator-valued observables, approximating their quantum 
counterparts in Hilbert-Schmidt topology. 

The proof makes use of a pointwise-positivity manipulation, which seems necessary in working with the 
norm and the precise form of the nonlinearity. We employ the Husimi function as a pivot between the classical 
probability density and the Wigner function, which - as it is well known - is not pointwise positive in general. 
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1. Introduction 

1.1. Formulations of the problem. The time evolution of a density matrix A'^{t) in a self- 
consistent field of Hartree type is described by the initial value problem 

ted,A%t) = [H^,A^{t)], 

A%0) = Al ^ ■ > 

where [-.■] denotes the commutator, e := h is the Planck constant and Aq is the initial datum. 
The dimension of the system is chosen equal to three. 

We will be interested in initial data of the mixed state form; see the statement of Theorem 
12. II for details. We are going to look at the dynamics (11.11) in the regime e <^ 1, the so-called 
semiclassical regime, where we expect the dynamics to "approach" the one of the corresponding 
classical system. The Hamiltonian for a particle of mass m is given by 

= A + V (1.2) 

and \^ is a self-consistent potential determined by 

V{x,t)= [ ^{x-x')pl{x')dx', (1.3) 



where 4>{x) is the two-body interaction and Pt{x) = p^{x,x,t) is the position density given 
in terms of the integral kernel p'^{x,y,t) of the density matrix A'^{t) (see below). The pair 
interaction potential </> is assumed to be spherically symmetric, a natural assumption from the 
physical point of view. 

It will be shown that the Wigner function corresponding to the operator (see below) 
converges in to the probability distribution g = g{x,k,t) that solves the corresponding 
classical Vlasov equation: 

dtg + k- d,g = T^g, (1.4) 

where 



{Tog) {x,k,t) = {d^(t)'k pI) ■ dkg{x,k,t) = yj ^^V 9x4>{x - y)pl{y)j ■ dkg{x,k,t), 

(1.5) 

and pf = pI{x) is the spatial probability density associated with g{x, k,t), namely: 

p?(x)= f dkg{x,k,t). (1.6) 



Moreover, a semi-classical approximation of A^ will be constructed; the precise statement of 
the result is given in Section 12. 1[ 
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Problem (11 .ip has several equivalent formulations, each of them yielding a corresponding 
initial value problem. For example, by looking at the integral kernel p'^{x,y,t) defined by 

{A%t) G){x) = j p'{x,y,t)G{y)dy, (for any G G L\R^)) (1.7) 

it is easy to check that it satisfies equation (13. 2 p - sometimes called von Neumann equation. 
Another description is in terms of the Wigner function f^{x, k, t) which is related to the density 
matrix A'^(t) through the following (Weyl) transformation: 

{A%t)G){x) = I f(^y,k,t)e^(-y)'lG{y)dkdy, (for any G e L\R')) (1.8) 

(see equations ( I1.12p . (13. 6 p for the relation between the Wigner function and the Weyl symbol). 
The Wigner function satisfies equation (13.90 - sometimes called Wigner equation. A very 
straightforward connection between the three descriptions is that 

3 

M^II/fS = ||P^(-, ■)I|l2(R3xR3) = (27re)2 ||/=||^2(lR3xM3), (1.9) 

where ||-||j|^5 is the Hilbert-Schmidt norm. This in particular allows to translate easily 
estimates between the different formulations, and to transfer approximations from the function 
level to the operator level. 

It must be noted that a natural assumption for the initial datum in the Density Matrix 
Formalism is that the operator is positive semi-definite and trace-class (with trace equal to 
one), i.e. it has a singular value decomposition of the form 

= ^ K\Um) {Um\ (1-10) 

with I Iwrral |l2(]k3) = 1, ^ 0, Xm = 1- It is well known that the trace is preserved in 
time. This is an important physical fact, as the trace is the quantum counterpart of the total 
probability of a classical density in a statistical formulation. The trace is given by tr{A'^) = 
J p'^{x,x)dx = J f^i^x, k)dxdk (see also Lemma \3A] and the discussion below). 

1.2. Physical context. Eq.n (II. ip describes the situation in which we have a large number of 
particles in a mean- field regime (see for instance [SH [ISl [S] , for the case of smooth potentials, 
and [H [121 El [HI [IHl [28] for more singular interactions). In this paper we want to study the 
semiclassical behavior (for e — 0) of the solution of equation (II. ip . 

The Wigner function is a well known tool in the study of the semiclassical limit of quantum 
dynamics (see definitions (13. 4p and (13. 5p below). Indeed, there are many works using the Wigner 
function to study the semiclassical limit of a number of problems (linear, non-linear, stochastic, 
systems etc) see e.g. [221 [201 [131 [23 [30] for a very small selection, and the references 
therein. One of the main advantages, is that the (formal, at this level) limit, as e — )■ 0, of the 
Wigner equation is typically some familiar equation of classical statistical mechanics. In that 
sense, the "correspondence principle" between classical and quantum mechanics is quantified in 
a straightforward, easy to present way. Indeed, for the problem we deal with here, for example, 
guessing the result from formal calculations is pretty straightforward. 
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In most of the existing literature, the notion of convergence is in weak topology (see e.g. 
the works mentioned above for precise statements). Indeed, the weak-* semiclassical limit for 
this problem is worked out e.g. in [20]. In fact (outside coherent states techniques), until very 
recently, virtually all the results were in weak topology. There is a natural analytical question 
of understanding when and why convergence in some natural strong topology fails; moreover, 
if one has possible numerical applications in mind, it would be desirable to know e.g. whether 
(possibly large in or pointwise sense) oscillations develop or not. Quantifying constructively 
the rate of convergence in terms of the data of the problem is also another natural question. 

Another big family of methods that yield strong topology semiclassical asymptotics (in 
linear as well as nonlinear problems) is based on coherent states (e.g. [HI [13 HSl [ISl 123 El E])- 
However, this is not really pertinent here, as coherent and mixed states (the kind of data we 
treat here) in nonlinear problems behave quite differently. 

In the quantum-classical correspondence, the idea is that the Wigner function converges, 
in the semiclassical limit, to a classical phase-space probability measure. However, the Wigner 
function itself is not pointwise positive in general. Working around this fact will be among the 
main points of the proofs. Indeed, it has often been remarked that the extensive arsenal of 
positivity techniques, developed in the context of classical phase-space equations, would be a 
good ingredient to transfer to their quantum counterparts. This becomes even more important 
in non-linear problems. In this paper we employ such pointwise-positivity techniques to the 
(nonlinear) Wigner equation, for the first time to the best of our knowledge. The key idea is 
to work with the Husimi function, a variation of the Wigner function, which does translate 
the operator positivity into pointwise positivity on the phase-space. The equation for the 
Husimi function itself has been derived only recently in closed form, and it is of infinite order 
in general [1]; it helps us in guessing the precise manipulations that are needed here. However, 
once formulated, the estimates we need can be proven without using the infinite-order Husimi 
equation itself. It should also be noted that the Husimi function is used only in the proof; 
it does not appear in the statement of the result, which is formulated in terms of the Wigner 
function. We believe that such use of positivity techniques in nonlinear Wigner equations could 
provide a fruitful approach in other problems as well. 

1.3. Notation. We specify that here and henceforth we use the following conventions for the 
Fourier transform in M"^: 



The Weyl Quantization is defined as follows: for any F e L'^iM? x M?) and G G L'^{M?), 
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We denote by H'^{W^) the Sobolev space of functions in L^(]R'^) whose derivatives 

up to the order u are also in L^(]R'^), i.e., for any function G on M'' 

V 

\\G\\h^I^9!^) = l|f^x^(^)lll,2(R<*) (1-13) 

|a|=0 



We denote by C^{R'^) the space of continuous and uniformly bounded functions whose 
derivatives, up to the order k, are also continuous and uniformly bounded. 

1.4. Organization of the paper. The plan of the paper is the following. Section 2 is devoted 
to the statement of our main result together with some remarks concerning it and fixing the 
notations. In Section 3 we recall the main features of the Hartree dynamics rephrased in the 
Wigner formalism. In Section 4 we introduce the definition and various properties of the Husimi 
Transform (HT). Then, in Section 5 we recall the main features of the Vlasov evolution, namely, 
the classical dynamics we recover in the limit e — 0. Finally, in Section 6 we prove the main 
result of this paper. 



2. The main result 

2.1. Statement of the result. All over the paper we make the following assumptions for the 
interaction potential that we suppose to be spherically symmetric: 

z) G H^R^) (2.1) 
ii) [ c/^ 10(^)1 < +00, n = 0,l,...,4. (2.2) 

Theorem 2.1. Under the above assumptions on (p, let A'^{t) be the solution of the Hartree 
problem and f^(t) the Wigner function associated with it. Denote f^i^O) = /q, and let us 

suppose that: 

• 3 C,C' > such that 

II/oIIh3(m3xk3) <C, J f^ix,k) \k\^dxdk < C and J f^{x,k)dxdk = 1. (2.3) 

• 3 Mq > 0, a G (0, 1] (independent of e), such that 

If^l'dxdk = 0(e^°), (2.4) 



\f^\dxdk = 0(e"), (2.5) 



Moreover 

||/o^||Li(R3) = 0(e"+ie*^o/i6.). (2.6) 
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Let {qq}^ be any family of probability distributions bounded in H^{M.^ x M^) (uniformly in 
e), supported on the set {{x,k) eM? ^M? : < Mq}, such that 

ll/o " fi'ollL2(]R3xiR3) = 0(5"), (2.7) 

(for example = XMofo / I XMofodxdk where /q is the Husimi Transform of fQ as in ( [^. j| j and 
Xmo is a smooth function identically equal to 1 for \k\ < Mq/2 and vanishing for \k\ > Mq, see 
Proposition \2.1\ below). 



Then, if we denote by g^{t) the solution of the Vlasov equation ljil-4\ ) with initial datum 
g^{0) = go , there exist positive constants Co,Ci,C2 (else we need to say independent of e, t), 
such that 

lirW -/WIL2(M3,M3) ^ C^e^^^^'^'eh. (2.8) 
In particular it follows that the density matrix A^{t) can be approximated by the semiclassical 
operator B"^ (t) whose Wigner function is g'^{t). More specifically 

M!W^^!Mk?^C„e---y" (2.9) 



HS 



In the statement of the theorem we do not require the existence of a semiclassical limit 
/o ~^ 9oy strictly need only assumptions on the quantum data /g. These assumptions 
guarantee the existence of semi-classical initial data g^ fulfilling the hypotheses of Theorem 
12. H as established by Proposition 12.11 below. On the other hand, if one supposes a priori the 
existence of such a family, the statement of Theorem 12.11 holds, without assuming conditions 
dUD, iUM and ([21D . 

Proposition 2.1. Let xmq = Xa/o(I^I) ^ monotone C°° function satisfying xmq = 1 ^/ 
\k\ < ^ and XMq = if \k\ > Mq, and let /q be the the Husimi Transform of fQ as in 
Suppose gl be defined by 

gl = ""'^f , (2.10) 
J XiviJodxdk 

then 

(1) J gQdxdk = 1 and g^ > 

(2) M\h-^=0{1) 

(3) suppg'QC{\k\<Mo} 

(4) \m-9'o\\L^=0{e-). 



The proof of Proposition 12.101 is given in the Appendix. 
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Remarks: 

• For a sharper expression on the behaviour in time of the error see (the end of) the proof 
in Section [61 In this form, the constants Cq, Ci depend on the i^^-norm of the initial 
data /q and Qq, the initial total energy of /q, suitable moments of (as in equation 
(12. 2p ). and All these quantities are bounded uniformly in e by assumption. The 
constant C2 depends on those quantities mentioned above that involve only 0. 

• The investigation of the semiclassical limit of the Wigner Transform by looking at the 

asymptotics arises quite naturally because such a norm is invariant under the time 
evolution. However, while for the linear case everything goes on easily (provided that 
the potential is sufficiently smooth; see e.g. [211 [29| [2]), for the nonlinear case one has 
to face an extra difficulty, which is the motivation of the present paper. In fact the 
L^-norm of the difference between the Wigner Transform and its classical counterpart 
is estimated in terms of the L^-norm of the same difference. Therefore, to conclude, we 
need a control of large momenta. This could be achieved by the energy conservation, but 
for an effective use of it we would need the positivity of the Wigner Transform, which 
is not the case. This difficulty has been overcome by using the Husimi Transform. 

• For the sake of concreteness we work in dimension three, but our results hold as well in 
any dimensions. We also expect the result to hold without any important differences if 
a regular enough external potential is added. 

2.2. Remarks on the Initial Data and Regularity Assumptions. 

Remark 2.1. The assumptions we made on cf) and guarantee the existence and uniqueness 
for the Vlasov equation ([i.^[ ) in the space of probability measures. 

Remark 2.2. An explicit example for which the assumptions of Theorem \2.1\ are verified, 
is a superposition of coherent states, namely, for an e— independent probability density go G 
H^{R^ X M^) and supported on {(x. A;) G x - |^| < Mq}, 

f^{x, k) = j dx'dk'Seix - x)5e{k - k')goix', k'), (2.11) 

where S^ix — x')6f.{k — k') is the Wigner transform of a coherent state centered in (x', k'). (see 
for instance [2Z]j- We observe that in this case the exponent a in \2. 7\) is equal to 1. 



In this case, any of the families {f^ole Theorem \2.1\ will converge to go in L^(M^ x M'^). 

Remark 2.3. Consider the case of a superposition of coherent states discussed in Remark lKE . 
For a pure coherent state centered in {x', k') and described by a density matrix A^,, ^, and a 
Wigner function f^'^k'i^^^)' have 

\Wx',k'\\HS = 1, il/x',fe'IU2(R3xR3) = — —j^. 



In contrast, when we deal with a mixture of coherent states (as in 1^2. we have ||/o||l2 < C 
(C independent of e). As a consequence (see U.9\) ). the corresponding density matrix has 
vanishing Hilbert- Schmidt norm. This is the reason why we consider the relative error in Ii2.9\) 
(the situation we take into account is the mixed state one). 

3. The Hartree dynamics in the Wigner picture 
The semi classical Hartree equation (with unit mass, m = 1) for a pure state is 

iedtu' = —u' u'. (3.1) 

For a mixed state, we have to pass to the von Neumann equation for the kernel p'^{x, y, t), 



+ (0^pf) (x) - (-f A, + (0^pf) (y))] p^ (3.2) 



where we put: 

pUz) ■.= p%z,z,t), (3.3) 

namely pl{z) is the spatial probability density associated with the quantum state described by 
p%x,y,t). 

The Wigner Transform of the wavefunction u^{x,t) is defined by, 

W'[u']{x,k,t):= Q-^ j e^^'^n^(a; + £|,t)n^(x-£|,t)d2/. (3.4) 

More generally for a mixed state described by a density matrix, namely a positive trace class 
operator A^(t) with kernel p'^{x,y,t), the Wigner function is 

3 r 

e'y^p'{x + e-,x-s-,t)dy. (3.5) 



/.(,,M) = (^) / 



The Wigner function is intimately related to the Weyl symbol, but one should be cautious with 
the scaling in the Planck constant: 



{A'{t) G){x) = j r (^^, fc, t)e'^'^-y^-^G{y)dkdy = Op^'y'{e^f')G{x) for G e 

(3.6) 

The spatial probability density can be easily expressed in terms of the Wigner function as 
well, namely 

pl{z):=p%z,z,t)= [ f%z,k,t)dk, \f zeR'. (3.7) 

Defining 

V{x,t)= I (t){x-x')pl{x')dx' = [ (pix-x') I f{x',k,t)dkdx', (3.8) 
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the self- consistent Hartree potential, the Wigner function evolves according to the Wigner 
equation, 

dtr + k-dj^ = Tfr, (3.9) 

where, for any w G L^(]R^ x R^) 

?■ r^/^ r 

iTiw) (x, k) = — — ^ / dA / d^ V{S, t) e'^^ {S ■ dk) w{x,k + e\S, t), (3.10) 

and V{S, t) is the Fourier transform (with respect to the space variable) of the potential V{x, t) 
defined in (13. 8p . It will be useful to observe that (13. 8p implies 

V{S,t) = ^{S)n{S). (3.11) 



It is well known (and easy to check), that the dynamics (13. 9 p preserves the integral of the 
Wigner function on the phase space x i.e. the trace (see [7] and Lemma [3.11 below) . 
This corresponds to the conservation of the L^-norm of the wave function u^, in case of a pure 
state, or to the conservation of the trace in case of a density matrix. For this reason we will 
have 

/ / f{x,k,t)dkdx= p'{x,x,t)dx = l, for any t > 0. (3.12) 

Lemma 3.1 {L^ regularity). Consider the initial value problem for equation liS. 9\) with initial 
datum /q. Under our assumptions on (p (more precisely, it is sufficient that (p G L^(R^)^, 
the trace is preserved by the time evolution. Moreover, if /q G L^(R^ x R^) , /^(t) stays in 
L^(R^ X R3) for allte R. 

Proof 

The trace associated with a Wigner function /^(t), namely, /(/^(t)) = tr (Op'^'^^'(£:^/^(t))) , is 
easily seen to be 

I{f^[t)) = j J f'{x,k,t)dkdx = j pl{x)dx (3.13) 

where the dx integral is understood to be absolutely convergent. The result for the preservation 
of the trace itself can be found in |^ (under the assumption (p G L°°, which obviously holds 
in this context). So far the dk integral does not have to be absolutely convergent, but only 
Cauchy-PV. 

Now, since pi G L^(R'^) for any t (and in particular ||pf ||2,i(ir3) = IIPoIIli(ir3) ~ -'-)' follows 
that pf G L°°(R'^) for any t (indeed ||pf ||/^oo(]83) < 1), and therefore 
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Now we can rewrite the Wigner equation (13.91) as 

{dt + k. d.)f = j^J i/(g)e-^^ ^^^^'^ + ^'^^-^^^^'^-^'^^ d^. (3.14) 
[Ziry J e 

One readily observes that the L^-norm (with respect to x and k) of the rhs of (I3.14p is bounded 
by 

2 

^^^;^||0||li(r3) ||/1Ui(r3xir3). (3.15) 

Since the free propagator (associated with the Ihs of fl3.14p ) preserves the L^-norm, the result 
follows by applying the Gronwall lemma. Observe that the constant grows like e~* (i.e. di- 
verging behaviour as e — 0) but this does not play any role here and we get all we need in 
justifying that the phase-space integral of is absolutely convergent. 

□ 



In the present context we need extra regularity properties in the framework of the Wigner 
formalism. Actually we can establish the following 

Lemma 3.2 (Sobolev regularity). Assume the potential cf) to satisfy the condition 

J dS\^{S)\ \Sr+' <+oo, 
for some m > 0. Then, for any T > 0, there is a constant C such that 

ll/^(^)IU™(R3xR3) ^ ||/^(0)||H™{R3xR3)e'^'^ (3.16) 

In particular for m = (i.e. looking at the L'^-norm), we have 

||/^(7')||l2(M3xM3) = ||/^(0)||l2(k3xR3). (3-17) 

Proof: By using the same observation as before, i.e. that |V^(5')| ^ ||l°°|(/'(5')| ^ |0('S')|, 
it follows that equation (13.91) can be treated in the same way as a problem with a smooth, 
time-dependent potential. The proof for the corresponding time-independent linear problem 
[29| [2] can be adapted to that end in a straightforward manner - this can be seen in more detail 
e.g. in [26] . 

4. HusiMi Transform and Husimi dynamics 

Given a Wigner function associated to a physical state (pure or mixed), we define the 
Husimi transform (HT) as 



r (x, k) = -F„-:,^.,,[e"i(-^+^''^) J-.,,^.,4/(., y)]] = 

= 7 — rr / e e E f[x,k)axak 
[Trey Jm3xr3 
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We sometimes denote by $ the smoothing map 

/ ^ $(/) = /. 

The remarkable feature of the HT rehes on the fact that f^>0 (see also the discussion 
in the Introduction). Indeed, in case of a pure state described by a wave function u, it can be 
verified by direct computation that 

[-Key J]R3xK3 

^ 0. (4.2) 



'XJl 

2 



{tie)* 
More generally, 

Je^^^k) = (4.3) 
(vre) 2 

where is the density matrix associated with the Wigner function 
Applying the map $ to the Wigner equation (13.91) . one finds [1] 

dtf + [k-d, + ^-d,- dk) r = ff r , (4.4) 

where 

?• /"^/^ r ~ F 

Tlw = —— d\ dSV{S,t)e'^^{S-dk)w{x + i-S,k + e\S). (4.5) 

The key observation is that, up to a small error, equation (14. 4p can be recasted as 

dtf + k-dJ' = Tjr + Eit), (4.6) 

where 

E{t) = E,{t) + E2{t) (4.7) 

and _ _ ~ _ ~ _ 

E, = -| 5x ■ 4f , i?2 = T// - T/r = <|.(T//) - T/r. (4.8) 

It is straightforward to observe that, according to (I3.10p (and f IS.lip ). we have 

yiw^ (x, k) = — ^ f ^ dX I dS V{S,t) e^^^ {S ■ dk) w{x, k + eXS,t), (4.9) 



1/2 



where 

V{S,t) = e-^mt) = HS)%{S) = HS)p{{S), (4.10) 
and the last equality follows easily by direct computation by setting: 

p{{x)= [ dkf%x,k,t). (4.11) 
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This observation will lead to the proof of Lemma 16.11 

The point here is that, up to a small error, the (non- negative) Husimi function satisfies the 
same nonlinear equation with self-consistent potential as the Wigner transform does. Indeed 
this observation was one of the main findings of [H [2], namely that (at least formally) to the 
leading order the Husimi equations are like the Wigner equations, but the potential has been 
replaced by a mollified version of itself. That is used here, since the dk marginal of the Husimi 
function is the mollification of the marginal of the Wigner function. Thus we preserve the 
structure of the quantum phase-space equation, while we change our function with one that 
remains non-negative. 

It must be noted that, once we formulate E2 = ^{Tl f) — T/ we do not really need the 
infinite order machinery to proceed (to the proof of Lemma [6.11 in this case). 



The Vlasov equation describes the situation in which we have a large number of classical 
particles in a mean- field regime (see for instance [H |23l |32] , for the case of smooth potentials, 
and [17], for more singular interactions). 

Denoting by $y(x, k) the flow associated with the system: 



one can easily verify that the solution g'^{t) of (11.41) with initial datum (see the claim of the 
Theorem 12.11 ) is obtained by propagating the initial datum through the characteristic curves 
of the flow $^(x, /c), namely. 



Therefore in proving existence and uniqueness of the solution of (11.41) one has to deal with a 
system of ODEs with a self-consistent field (see (15.11) ) and the smoothness of the potential is 
sufficient to apply a fixed point argument (see [H [10], [2l] ) . 

Also for the Vlasov equation we need some regularity properties of the solution and we 
will make use of the following 

Lemma 5.1 (Vlasov regularity). Assume the potential G C'^'^^iE?) for some m >1. Then 
there is a constant C such that 



5. Vlasov evolution 




(5.1) 



g'{x,k-t)=gl{^y' (x, k)). 



(5.2) 



llfi'^(^)llH'"(M3xIR3) ^ ||fl'oll/^'"(IR3xR3)e * 

The same proof holding for the Wigner case does apply here. 



(5.3) 



6. Proof of Theorem 12.11 



In the course of the proof we will denote by C any positive constant, possibly depending 
on 0, /q or 5fQ, but neither on t nor on e. 



13 



It is not difficult to show that under suitable smoothness assumptions the Wigner and 
Husimi functions are close in L?. In particular, by direct computation the following inequality 
can be proven 



nt) - nt) 



L2(R3xI 



(6.1) 



(see also Lemma A.l in [2]). 



Moreover we can also show that practically solves the nonlinear Wigner equation up to 
a small error in (see fl4.6p ). This will be used in the main body of the proof below. 
Indeed, we recast equation (14. 6 p as 



dtr + k-dj'-T!r = E{t), 



(6.2) 



(see Section H]) and observe that (16. 2p can be seen as a Wigner equation with a time-dependent 
source term. The error E{t) is ensured to be small by the following: 

Lemma 6.1. Assume /q to satisfy all the assumptions of Theorem \2.1\ and E{t) he defined by 
equations ( [^.6| j-( [7^ . Then 

<Ce^'e. (6.3) 



\\m\ 



Proof 

We ffist bound Ei (see (gTI) and (jM])). We have 



(6.4) 



where we estimated, uniformly in e, the H'^-norm of with the if^-norm of and we used 
property (13.161) for 

Moreover, since we have 



E2 = Tif-Tif 



(6.5) 



then, by (I4l9|) and (143!) we get 

•1/2 



(27r) 



1/2 



dSV{S,t)e'^''{S -dk) \J%x + i-S,k + e\S)- f%x,k + e\S)^ (6.6) 



and we remind that E2 = E2{x, k,t). Thus, by taking the Fourier transform J^x,k^p,q we find: 
.1/2 

dX I dS (J)(S)p{(S)e'''''"' S ■ (ia) fHv- S.q) et^'^e-t^*' - 1 , (6.7) 



namely 
E2{p,q) = 



(2vr)= 



1/2 



1/2 



dX 



1/2 



dS4){S) p{ (S) e''^^'^ S ■ (iq) f{p - S, q)e-iP' 



1 — e 2 e2 
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By applying the Taylor formula, for some C, G (0, e) we get 



and hence 
Finally we obtain 



F 9 F 9 



1 - e ^s'^e^P^ 



<-{S' + \p\\S\)<e{S' + \p-S\\S\). 



\E,ip,q)\' < Ce' \q\' \S\ iS' +\p- S\\S\) Ifip - S,q)\ 

where we used f l6.10p and the uniform L°° control on p{ . 
Then, by applying the Cauchy-Schwarz inequality 



(6.9) 



(6.10) 



(6.11) 



E,{p,q)\'<Ce'\q\U / d5|0(5)||5 



dS\<P{S)\{S' + \p-S\\S\r\r{p-S,q)\' 

(6.12) 



Therefore 



(6.13) 



where we made use of property fl3.16p for and the assumptions we did on (p (see (12.21) ). 

Then, by (16.41) and (16.131) we obtain that (16. 3p holds true and the proof of Lemma 16.11 is 
concluded. 

□ 



Before moving on, let us make a remark on notation: as it has been explained previously 
(see e.g. the statement of Theorem 12. H or the discussion in Section 12. 2p the classical initial 
data may or not depend on e. However all the characteristics of interest here (Sobolev 
norms, support etc) are bounded uniformly in e. Consistently with that, and for simplicity in 
the notation, we will drop the superscript e, and refer to the initial data go and the classical 
solution g{t) in the sequel. 

For convenience of the reader, we recall that, according to the statement of Theorem 12. H 
5^0 satisfies the property: 

supp go C {(x, fc) G X : \k\ < Mq}, (6.14) 



for a certain constant Mq not depending on e. 
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By virtue of (16. ip . in order to prove Theorem 12.11 (more precisely, formula (12. 8p ) it suffices 
to bound in the remainder ht{x, k) = /^(x, k, t) — g{x, k, t); indeed we will show that 

||/it||L2(R3xR3) = writ) - (7(t)||L2(K3xM3) ^ Ce^^''*£^". (6.15) 

Then, (12. 8 P will follow straightforward by joining (16. ip and (I6.15p . 



(6.16) 



Proof of equation ( I6.15P By (ll.4p and (14. 6 p the evolution of ht is given by 

dtht + k ■ d,ht = Tjht + (t/ -Tl^g+ (tJ - To^) g + E{t). 
ho{,x, k) = f^{x, k) - goix, k). 
Let be the (L^ preserving) Wigner-Liouville flow associated with the equation: 

dth +k-d,ht=Tjht. (6.17) 

Then, from (I6.16P we have: 

ht = n'ho+ f ds n'~'{E{s)) + [ ds n'^' {ri{s) + r2{s))) , (6.18) 
Jo Jo 

where (see flXTOj) and ffT3|) ) 
ri(s) := (tI -Tpj g{x,k,s) = 

= — — / dA / dS<P{S)pi{S)e''- {S-d,)[g{x,k + e\S,s)-g{x,k,s)], 

(6.19) 

and 

:^ (tI _ rs) .) . (l d, a.,i. - ,) [p[(,) - p,(,)] ) . 

(6.20) 

Next, from (16.180 . we have: 

II^4l2(R3xR3) < ||^o|lL2m3xR3) + / ds \\ E {s)\\ ^2 m3 y^^3\ + 

Jo 

+ ^ ds (||ri(s)||^2(K3xR3) + lk2(s)||i2(K3xR3)) , (6.21) 

and, by Lemma [6. II 

L2(R3xR3) < II^0|Il2(R3xR3) + C {^''^ ^ ^ ^ ^'^^ ( 1 1 ''l (■^) 1 1 (R3 xR3 ) + 1 1^2 (s) 1 1 ^2 (ir3 xR3) j • 

(6.22) 



16 



Moreover by (12. 7p and estimate (16.11) . we easily get: 



fo - 90 



< 



L2(-k3xI1 



fo /( 



L2(IR3xIl 



+ II /o 5'0|li2(]u3xM3) 



(6.23) 



Then, we finally get: 



* IIl2(r3 



< C(e^* + 1)5" + l^ds (||ri(s)||^.(j,3xM3) + ||r2(s)IL.(M3xM3)) • (6-24) 
Next we evaluate the L^-norm of ri(s). 

We observe that, by virtue of the positivity of the L°°-norm of p{ is uniformly bounded. 



In fact, we have: 



< 



Pt 



dx 



dk f{x,k,t) 



(6.25) 



where the last equality is easily obtained by direct computation (see (14. ip and (I3.12p ). Then, 
by applying the Taylor formula in (I6.19p . we can estimate the L^-norm of ri(s) as follows: 



<Ce' 



dS\<t>iS)\\S\' ||^7(.)||^3(M3xM3)<C^ 



(6.26) 



where, in the last inequality, we used property (15.31) for g{s) and assumption (12.21) for 0. 
Then, by ^7M^ we get 

ft 



\h 



t\\Ly 



<C(e^* + l)£" + C(e^* + l)£2+ / ds \\r2{s)\\^,, 

Jo 



(6.27) 



Now, let us look at the L^-norm of r2(s). 
By (K^ we have 



\r2{S)\\L2i 



< 



< dx dk \dkg{x, k, s) 
We split the integral 



dy \d^(f){x-y)\ / dw f'{y,w,s)-g{y,w,s) 



dw 



fs{y,w) -gs{y,w] 



(6.28) 
(6.29) 



into the two domains \w\ < M and \w\ > M, where M is chosen in the following way. 
If X{t), K{t) is the classical flow generated by the force field 

dyd^(t){x - y)Ciy,t), 
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where C is any spatial probability density, then 

\K{t)\<\K{0)\ + \\dML^t. 



Therefore, by virtue of assumption ( 16.14^ on the initial datum qq, there exists a positive constant 
M for which (for t < T arbitrary but fixed) 

g{x,k,t) = if |A;|>M. (6.30) 

Clearly M = M{t) depends on time and it is straightforward to check that: 



M{t) = Mo + 



(6.31) 



where Mq is the same as in fl6.14p . 
Hence, by fl6.28p . we have: 



lk2(-S)||^2(K3xR3) 



< dx dk \dkg{x, k, s)\ 
dx / dk \dkg{x, k, s)\ 
+2 / dx / dk\dkg{x,k,s] 
X (^3^^ \dx4>{x-y)\ J 
2 dx dk\dkg{x,k,s 



\w\<AI 



w\>M 



\w\<M 



dw 



|tu|>A/ 
\ i2 



-2 / dx dk \dkg{x, k, s)\' 



dy \d^(f){x-y)\ 
dy \d-,(j){x - y)\ 
dy \d^(f){x - y)\ 
f{y,w,s) - g{y,w,s 
dy \d,^(f){x-y 
dy \d^(f){x-y)\ 



dw 



dw 



dw 



f{y,w,s) - g{y,w,s] 
r{y,w,s) - g{y,w,s) 
f{y,w,s) - g{y,w,s) 



+ 



+ 



< 



dw 



\w\<AI 



dw 



'\w\>M 



f%y,w,s) - g{y,w,s 
f{y^w,s) - g{y,w,s) 



(6.32) 



We first bound the first term on the right hand side of (16.321) . We find 



dy \d^(l){x-y)\ 



dw 



w\<M 



fiy,w,s) - g{y,w,s] 



dy / dw \d^(f){x-y)\ X\w\<m{w) f{y,w,s)-g{y,w,s) 



< 



r{s)-g{s) 



..(M3.M3)^r^'ll^-'^"^^ II^Ji^(M3xM3). (6.33) 
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Therefore, we obtain 



2 dx dk\dkgix,k,s)n / dy 



dw 



< Ce^' M'^ \\h 



Ce^%C + sf\\h,\\' 



\w\<M 



f'{y,w,s) - g{y,w,s) 



< 
(6.34) 



by virtue of Lemma 15.11 ( 16.31^ and assumption (12. ip on 0. 
Now let us look at 



2 dx dk\dkg{x,k,s)f ( / dy \d:,(j){x - y)\ 

By fl6.30p we have 
2 J dx dk\dkg{x, k, s)\ 

2 J dx dk\dkg{x,k,s)\ 

< A J dx dk\dkg{x,k,s)f' 

+4 / dx dk\dkg{x,k,s)f 



dw 



|tu|>A/ 



f{y,w,s) - g{y,w,s] 



(6.35) 



dy \d^(j){x - y)\ 
dy \d^(j){x - y)\ 
dy \d^(f){x-y)\ 
dy \d^(j){x-y)\ 



dw 



w\>M 



f%y,w,s) - g{y,w,s] 



dw 



w\>M 



|ui|>e 



f%y,w,s] 



< 



dw 



dw 



M<\w\<e- 



f'{y,w,s) - g%y,w,s] 



(6.36) 



where a; > will be fixed later (see f l6.49p )). Here g^{y, w, s) denotes the solution of the classical 
Liouville equation generated by the force field 



namely 



dyd^(t){x - y)p^{y,t), 



(6.37) 



and the initial datum is the same of the Vlasov evolution, namely, g^{y., w, 0) = go{y, w). Note 
that g'^ enters freely in the game because it satisfies the support property (16.301) . 



Let us estimate the term: 



4 / dx dk\dkg{x,k,s)f ( / dy \d^(f){x - y)\ 



dw 



\w\>e- 



f%y,w,s] 



(6.38) 
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Note that here the positivity of is crucial because it allows us to use the energy conservation. 
Indeed 



<^\\9{s)\\li\\d,m^( [ dy I dwf{y,w,s)] < 

\Jm3 J\w\>e-'^ / 



< 



\w\>e~ 
2 \ 2 



<Ce'-[ dy dw — f{y,w,s)] < Ce^^ (6.39) 

To show the last inequality, we denote by ■y{w, s) = f dy /^(y, w, s) the distribution of momenta 
(which is obviously positive and with integral in dw equal to one) and by 7(w, s) the smoothed 
version of 7(w, s). Then, denoting by $' the smoothing acting only on the momentum variable, 
we set 7 = $'(7) (clearly, the action is exactly as in (14. ip for $). As we observed for the spatial 
distribution pi (see (14.101) ). even in this case it is straightforward to see that the smoothing 
commutes with the partial integration on the phase-space, namely 



7(u;,s) = $'(7)(u;,s) = $'(y dyf{y,w,s)) = j dy ^{f){y,w, s) = / dyf{y,w,s). 

Therefore, we can write the term J^^ dy J^^ dw ^ f^{y, w, s) in (16.391) as 

f f w'^ ~ f w"^ ^ 

dy dw — f{y,w,s) = I dw — -f{w,s). 

Jm.3 Jk3 Z J^3 z 



Moreover, we have 



I dw — 7(w, s) = / dw — e^'^7(w,s), (6.40) 
Jrs 2 J 2 

because the action of the smoothing operator is exactly the same of the heat flow e^^ (see 
(14.11) ). Now we use the well-known (and easy to check) property 



dw — e'^^'y{w, s) < / dw ( ^ ) ^) ~ / ~2~ ~^ 



to conclude that: 



w'^ ~ f w"^ f f w'^ 

dy I dw — f{y,w,s)< dw — 'y{w,s) + Ce= dy dw — f{y,w, s) + Ce. 

Z Z ./ttpS /ma / 



(6.41) 
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Now, by the energy conservation: 

[ dy [ dw%- f{y,w,s) < 

Jr3 Jr3 / 



< dy dwf{y,w,s) 



[ dy [ dwf^{y,w) 

JR3 Jr3 \ ^ 



* p-^ (0) 



c 



(6.42) 



where, in the first inequahty, we used that the potential is bounded from below ( as a 
consequence of assumption (12. 2p ) and the last bound follows from our assumptions . 

Thus, by f lCTj) and (l6^ we get ( ICTj) (the control on 9^0 is guaranteed by 
and the control on g is guaranteed by (15.31) ). 

Next we estimate the term on the last line of (16.361) . namely: 

2 



A dx dk \dkg{x, k, s) 



dy \dx(f){x -y)\ I dw 

'M<\w\<e- 



< 



Cs ^—3ui 



f'-t 



LH 



(6.43) 



To control the L^-norm of /^(s) — g'^i^s) we set pi = /^(s) — g'^{s). Then, by (14. 6 p and 
(16.37p . the equation for is: 



(6.44) 



dsPt + k ■ d^pl = Tfpl + [Tf - TfJ r{s) + E{s). 
Poi^^ k) = foix, k) - go{x, k). 
We proceed as before using that the flow generated by —k-dx + T/ is isometric in L^. Therefore 

IIp^IIl^ < WpIWl^ + dr\\E{T)\\L2 + J^' \\ {Ti-Tl^rir)\\L2dT (6.45) 

The first two terms on the right hand side of (I6.45p have been estimated previously (see (16.231) 
and Lemma 16.11) and they give rise to 



L2 < Ce" and 



dT\\E{r)\\L2 < Cie^'' + l)e 



(6.46) 



Moreover, the last term in (16.450 can be estimated exactly as the term ri(s) (see (16.191) and 
(I6.26P ) because the Liouville dynamics for g^ controls the if^-norm (that is finite and uniformly 
bounded in e at time r = since ^^(0) = go). Therefore we find: 



Thus, finally we get: 



To^-T/jr(r)|U^rfr<C(e^^ + l)e2. 
Ib^lU. <C(e^^ + l)5". 



(6.47) 
(6.48) 
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Therefore, by setting 



a 



(6.49) 



we find 

4 j dx dk\dkg{x,k,s)\^ [ / dy \d^(l){x - y)\ 



dw 



M<\w\<e'- 



f%y,u!,s) - g'{y,w,s) 



< 
(6.50) 



with —3a; + a > 0. 

In the end, by collecting (16.341) . (I6.39P and (I6.50p . we obtain that: 



, (6.51) 



that gives an optimal bound for w = y, namely: 



'■s|Il2(R3xR3) 



(6.52) 



and 



ds \\r2{s)\ 



<C (e 



(6.53) 



By (I6.27P and (I6.53p . we can finally control the L^-norm of /if, obtaining: 

I|/^*IIl2(m3xm3) < C ((e^* + 1)5° + (e^* + 1 + t)e'r'^ + ds {C + sf' e^^ 1^.(^3 xm3)) ■ 

(6.54) 

The, by applying the Gronwall lemma, we find that 

||/..|L.(m3xm3) < Cie^' + 1 + t).f°e/o^^ (c+.)3/2e-^ 

and hence (16.151) holds true. 

Finally, by virtue of (16. ip and (16.551) . we deduce (12. 8p and hence the proof of Theorem 12.11 
is concluded. 

□ 

Remark 6.1. We observe that, for example, in the case in which the initial datum is given by 

2 

/{2.11]) . the rate of convergence is estimated by 57. 
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Appendix A. Existence of semiclassical initial data: proof of Proposition 12.11 

Let g'o be defined by fl2A0D . Set N = N{e) := (/ f^{x, k)xMMk\)dxdk)-\ 
Property ([1]) follows by the positivity of the Husimi function and ([3]) follows by construction. 
Since 



where C depends on the L°° norm of xmo ^^id its derivatives, ([2]) follows as soon as we show 
that N = N{e) = 0(1) (see below). 

So it is only left to check (jl]). 

Observe that 

^ ||(l-iV;^,,J/„-||^. + ||iVxMo(/o^-7^)IU^ ^ (A.2) 

^\1-N\ ll/olU. + iV||(l - Xmo)/o^I|l^ + iV||/o^ - Ml^ 

Obviously, the estimate for the second term of the last line of flA.2p comes from assumption 
(12.41) . and for the third term from the fact that ||/q — /o||l2 = 0(£:||/o H/^a) as shown in lemma 
A.l of P]. To see that the first term is 0{e°'), i.e. that |1 — = 0{e"'), one has to observe 

that (denote for brevity 5^ (A;) = — ^-e ~) 

(7re)2 

|1 - = I / (7^XMo - mdxdk\ ^ 11(1 - XMomilL^ ^ 

^ / 6'{k-k')f^{x',k')dx'dk'dk= J 6'{k-k')f^{x',k')dx'dk'dk+ 

\k\>r 
Ifc'l > 



\k\>^ \k\>^ 



+ / S'{k-k')f^{x',k')dx'dk'dk^ J S'{k)dk J \f^{x',k')\dx'dk'+ (A.3) 



Ifc'l < ^ 



+ sup \S%k-k')\ ll/o^lU, =0(0 



lfcl>^ 

Ifc'l < ^ 



where we made use of assumptions (12.51) and (12. 6p in the final step. 
Now obviously, 

|l-i|^Oe" |iV- 1| ^ Oe°|iV| ^ 
|A^- IK 0£"(1+ |iV- 1|) \N - 1\ ^ = 0{e' 



(A.4) 
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